Environmentally-Induced Rabi Oscillations and Decoherence in Phase Qubits 
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We study decoherence effects in a dc SQUID phase qubit caused by an isolation circuit with a 
resonant frequency. The coupling between the SQUID phase qubit and its environment is modeled 
via the Caldeira-Leggett formulation of quantum dissipation/coherence, where the spectral density 
of the environment is related to the admittance of the isolation circuit. When the frequency of the 
qubit is at least two times larger than the resonance frequency of the isolation circuit, we find that 
the decoherence time of the qubit is two orders of magnitude larger than the typical ohmic regime, 
where the frequency of the qubit is much smaller than the resonance frequency of the isolation 
circuit. Lastly, we show that when the qubit frequency is on resonance with the isolation circuit, 
an oscillatory non-Markovian decay emerges, as the dc SQUID phase qubit and its environment 
self-generate Rabi oscillations of characteristic time scales shorter than the decoherence time. 

PACS numbers: 74.50.+r, 85.25.Dq, 03.67.Lx 



The theoretical possibihty of using quantum mechan- 
ics to manipulate information efficiently [l[ has lead, 
through advances in technology, to the plausibility of 
building a quantum computer using two-level systems, 
also called quantum bits or qubits. Several schemes 
have been proposed as attempts to manipulate qubits in 
atomic, molecular and optical physics (AMO), and con- 
densed matter physics (CMP). In AMO the most promis- 
ing schemes are trapped ion systems and ultracold 
atoms in optical lattices Q . 

On the CMP side, the pursuit of solid state qubits 
has been most promising in spin systems 0, @| and su- 
perconducting devices [g, 0, [1] ■ While the manipulation 
of qubits in AMO has relied on the existence of qubits 
in a lattice of ions or ultra-cold atoms and the use of 
lasers, the manipulation of qubits in CMP has relied on 
the NMR techniques (spin qubits) and the Josephson ef- 
fect (superconducting qubits). Integrating qubits into a 
full quantum computer requires a deeper understanding 
of decoherence effects in a single qubit and how different 
qubits couple. 

In AMO systems Rabi oscillations in single qubits 
have been observed over time scales of miliseconds since 
each qubit can be made quite isolated from its environ- 
ment [2], however it has been quite difficult to imple- 
ment multi-qubit states as the coupling between differ- 
ent qubits is not yet fully controllable. On the other 
hand, in superconducting qubits Rabi oscillations have 
been observed [6|] over shorter time scales (500ns), since 
these qubits are coupled to many environmental degreees 
of freedom, and thus require very careful circuit design. 
Furthermore, extensions to the multi qubit regime also 
require further integrated circuit designs. 

In this manuscript, we analyze decoherence effects in 
a single superconducting phase qubit coupled to iso- 
lation circuits. (Phase qubits Q are superconducting 
qubits dominated by the Josephson effect, in contrast to 
charge qubits |^], which are dominated by quantization 



of charge.) The coupling of the qubit to the isolation 
circuit is emulated by a spectral density with an intrin- 
sic resonance within the Caldeira-Leggett formulation of 
quantum dissipation. 

The circuit used to describe intrinsic decoherence 
and self-induced Rabi oscillations in phase qubits is 
shown in Fig. [1] which correponds to an asymmetric dc 
SQUID The circuit elements inside the dashed box 
form an isolation network which serves two purposes: a) 
it prevents current noise from reaching the qubit junc- 
tion; b) it is used as a measurement tool. 

The classical equation of motion for such a circuit is 

CoJ+— leasing -—I + / dt'Y{t-t')^{t') = (1) 

*0 *0 Jo 

where Icq is the critical current of Josephson junction J 
in FiglU and $0 = h/2e is the flux quantum. The last 
term of Eq. [1] can be written as iu!Y{u;)j{uj) in Fourier 
space. The admittance function Y{uj) can be modeled 
as two additive terms Y(lo) = Yisoi^j-') + Yinti^^)- The 
first contribution Yiso{'-o) is the admittance that results 
when a transmission line of characteristic impedance R 
is attached to the isolation junction (here represented 
by a capacitance C and a Josephson inductance L) and 
an isolation inductance Li. Thus, Kiso(w) = Z~l{u)) 

where Ziso{(^) = (iujLi) + [R^^ + iujC + {iujL)^^'\ ^ is 
the impedance of the isolation network shown in Fig. [T] 
The replacement of the isolation junction by an LC cir- 
cuit is justified because under standard operating con- 
ditions the external fiux varies to cancel the current 
fiowing through the isolation junction making it zero bi- 
ased Thus, the isolation junction behaves as a har- 
monic oscillator with inductance L which is chosen to be 
much smaller than Li. The second contribution Yi„t{^^) 
is an internal admittance representing the local environ- 
ment of the qubit junction, such as defects in the oxide 
barrier, quasiparticle tunneling, or the substrate, and can 
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be modeled by li„f (w) — (Ro+iujLo) ^, where i?o is the 
resistance and Lq is the inductance of the qubit as shown 
in Fig. m 



broadening controlled by T. In addition, notice that the 
dimensionless parameter F/il^ — LLi/ (2i?(ii + L)) « 
R/L is independent of C. Therefore, when there is no 
capacitor (C ^0), the resonance disappears and 
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reduces to a Drude term with characteristic frequency 
il2/2r Ki R/L. The internal spectral density of the qubit 



FIG. 1: Schematic drawing of the phase qubit with an RLC 
isolation circuit. 

Next, we use the Caldeira-Leggett formalism to de- 
scribe the coupling of the phase qubit to the isolation 
network, via the spin-boson Hamiltonian 0] 



H = 



E 



SB, 



(2) 



written in terms of Pauli matrices ai (with i = x,y, z) 
and boson operators bk and &^ . The first term in Eq. ^ 
represents a two-level approximation for the phase qubit 
(system) described by states |0) and |I) with energy dif- 

where Er ~ 



ference Tiuj^ — y^8EcEj{l 
(2e)2/2C is the capacitive energy, and Ej ~ (<i>o/27r)/co 
is the Josephson energy. The second term corresponds 
to the isolation network (bath) represented by a bath of 
bosons, where bk and are the annihilation and creation 
operator of the fc-th bath mode with frequency ujk- The 
third term is the system-bath (SB) Hamiltonian 



HsB = -a,n(l|7|0)^Afc 

k 



('t 



(3) 



corresponding to the coupling between the isolation net- 
work and the phase qubit which appears as dt' S(t — 
t')'^{t') in the classical equation of motion Eq. 

The spectral density of the bath modes J(a;) = 
A^^ (w — ujk) has dimensions of energy and can be 
written as J{lo) — LoYieY {lo){^q / 2tiY , which leads to the 
compact form J(w) = Jiso{^) + Jinti^)- The spectral 
density of the isolation network is 



Jiso ('^) = ( ^ 



(4) 



where a = l?l{{L L^fR) « [L/LiY/R is the lead- 
ing order term in the low frequency ohmic regime, O = 
yJ{L + Li) / {LLiC) ~ I/VLC is essentially the reso- 
nance frequency, and T = l/(2Ci?) plays the role of res- 
onance width. Here, we used Li ^ L corresponding to 
the relevant experimental regime. Notice that Jiso{uj) has 
Ohmic behavior at low frequencies limi^^o Jisoi^)/'-^ — 
(<i>o/27r)2(L/Li)2/i?, but has a peak at frequency with 
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l + co^Ll/Rl 



(6) 



is a Drude term with characteristic frequency Rq/Lq. 
Notice that Jint{uj) also has Ohmic behavior at low fre- 
quencies lim^^^o •^mt('jj)/^ — ($o/27r)2/i?o- In order to 
obtain the relaxation Ti and decoherence times T2, we 
write the Bloch-Redfield equations 



^nmPnrn 



%klPkl 



(7) 



kl 



for the density matrix pnm of the spin-boson Hamiltonian 
in Eq. ([2|) and ([3]) derived in the Born-Markov limit. Here 
all indices take the values and 1 corresponding to the 
ground and excited states of the qubit, respectively, while 
^nm = {En — Em)/^ is the frequency difference between 
states n and m. The Redfield rate tensor is 



J^nmkl — -L ijnnk ^ Irnnk ^ Irrm ' nrrfc' 

where repeated indices indicate summation, and 



(8) 



^ Imnk ~ 



dte- '^^ * {HsBAm (t)HsB,nk (0)) , (9) 



where j = 1,2 and 771 — Unk amd 772 — ^im- Under these 
conditions, the relaxation rate 1/Ti — ^ Rnnnn becomes 
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Muj, 



^j(^„,)coth(^ 
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(10) 



where M = ($o/27r) Co has dimensions of mass x area 
(or energy x time squared) and is refered to as the mass 
of the phase qubit with capacitance Co, while wqi is the 
qubit frequency. The interpretation of T-^^ is as follows: 
for the system to make a transition it needs to exchange 
energy E — Tiloqi with the environment using a single 
boson. The factor coth(?iwoi/^BT) = n{ujQi) + 1 + n{LLjQi) 
captures the sum of the rates for emission (proportional 
to n(woi) + 1) and absorption (proportional to n(u;oi) of 
a boson), where n(woi) = [exp(?lti)oi/fcsT) — 1] is the 
Bose function. 

In Fig. [21 Ti is plotted versus qubit frequency luqi for 
spectral densities describing an RLC (Eq. 2]) or Drude 
(Eq. [5]) isolation network at fixed temperatures T — 
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FIG. 2: Ti (in seconds) as a function of qubit frequency ujoi. 
The solid (red) curves describes an RLC isolation network 
with parameters R = 50 ohms, Li = 3.9nH[, L — 2.25pH, 
C — 2.22pF, and qubit parameters Co — 4.44pF, Ro — oo 
and Lo — 0. The dashed curves correspond to an RL isolation 
network with the same parameters, except that C = 0. Main 
figure (T = 0), inset (T = 50mK) with Q = 141 GHz. 
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FIG. 3: Ti (in nanoseconds) as a function of qubit fre- 
quency uoi. The solid (red) curves describes an RLC iso- 
lation network with parameters -R = 50 ohms, Li = 3.9nH, 
L = 2.25pII, C = 2.22pF, and qubit parameters Co = 4.44pF, 
-Ro = 5000 ohms and Lo = 0. The dashed curves correspond 
to an RL isolation network with the same parameters, except 
that C = 0. Main figure (r = 0), inset (T = 50mK) with 
Q. = 141 GHz. 



(main figure) and T = 50mK (inset), for Jintii-^) = 
corresponding to _Ro ^ oo. In the limit of low temper- 
atures [ksT /hujoi <^ 1), the relaxation time becomes 
Ti(a;oi) = Mujoi/J{i^oi)- From Fig. [2] (main plot) sev- 
eral important points can be extracted. First, in the 
low frequency regime (wqi ^ ^) the RL (Drude) and 
RLC environments produce essentially the same relax- 
ation time Timc{0) = Ti^rUO) = T^fl ~ [L^/LfRCo, 
because both systems are ohmic. Second, near reso- 
nance {ijjQi ~ n), Ti^nLc is substantially reduced be- 
cause the qubit is resonantly coupled to its environ- 
ment producing a distinct non-ohmic behavior. Third, 
for {loqi > il), Ti grows very rapidly in the RLC case. 
Notice that for uqi > the RLC relaxation time 

Ti^RLC is always larger than Ti^hl- Furthermore, in the 
limit of LjQi 3> max{f2, 2r}, Ti^rlc grows with the fourth 
power of ujQi behaving as Ti^rlc ~ 2^1,0^^01/^^^, while for 
ciJoi S> ri'^/2r, Ti^nL grows only with second power of 
LUoi behaving as Ti,rl ~ '^Ti^o^^UQi/ft'^. Thus, Ti,rlc 
is always much larger than Tij^l for sufficiently large 
luqi. Notice, however, that for parameters in the experi- 
mental range such as those used in Fig [21 Ti.flic is two 
orders of magnitude larger than Ti jil, indicating a clear 
advantage of the RLC environment shown in Fig [1] over 
the standard ohmic RL environment. Thermal effects 
are illustrated in the inset of Fig. [2] where T — 50mK is 
a characteristic temperature where experiments are per- 
formed . The typical values of Ti at low frequencies 
vary from IQ-'^s at T = to lO^S at T = 50mK, while 
the high frequency values remain essentially unchanged 
as the thermal effects are not important for Tiwoi ^ ksT. 

In the preceeding analysis we neglected the effect of 
the local environment by setting Yinti^^) — 0. As a re- 
sult, the low- frequency value of Ti is substantially larger 
than obtained in experiment 13, lH- By modeling the 



local environment with Rq = 5000 ohms and = we 
obtain the Ti versus luqi plot shown in Fig. [31 Notice 
that this value of Rq brings Ti to values close to 20ns 
at r = 0. The message to extract from Figs. [H and [31 
is that increasing Rq as much as possible and increasing 
the qubit frequency ujqi from G.lf2 to 2f2 at fixed low 
temperature can produce a large increase in Ti. 

Although the Bloch-Redfield equations described in 
Eq. ([7]) capture the long time behavior of the density 
matrix, they can not describe the short time behavior of 
the system in particular near resonance where lvqi « Ct, 
where the RLC spectral density is very large. In this 
case, only the environmental modes with ujk = ^ couple 
strongly to the two-level system, like a two-level atom 
coupled to an electromagnetic field cavity mode that has 
a finite lifetime. This is best seen by restricting the 
Hamiltonian described in Eqs. ^ and ^ only to bo- 
son modes with ujk ~ ^ ~ cuoi- In this case, it is best to 
rewrite the spectral density as 



<T=±1 UJ^ 



n + ir 



(11) 



where a has the same definition as in Eq. (|4]), and 
— n — r'^/il. This reveals a resonance a.t u = fl with 
linewidth T, such that Jint{^ — ^) can be neglected for 
any non-zero value of Rq, and J{oj) ~ Jiso{oj)- 

When ijjk ~ ^ ^ wqi, the Hamiltonian in Eqs. (0) 
and (O can be solved in the rotating wave approxima- 
tion using the complete basis set of system-bath prod- 
uct states iVo) = |0)s » |0)b; IV'i) = |l>s ® |0)b; 
iV'fe) = |0)s ® \k)B, where |0)s and |l)s are the states 
of the qubit and |fc)B are the states of the bath. Hence, 
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the state of the total system at any time is 

(f>{t) = coV'o + ci{t)'ilJi + ^ Ck{t)'ipk, 



(12) 



with probabihty amphtudes cq, ci(i), and Ck{t). The 
amplitude cq is constant, while the amplitudes ci{t) and 
Cfe(i) are time dependent. Assuming that there are no 
excited bath modes at t = 0, we impose the initial con- 
dition Cfc(O) = 0, and use the normalization — 1 
to obtain the closed integro-differential equation 



CiW = - / dh 
Jo 



f{t-h)ciih), 



(13) 



where the kernel is the correlation function 

huj 



fir) = / dioJiij) 



coth 



2kBT 



cos{ujt) + i sin(ti;T) 



directly related to the spectral density J{oj). 
present case the reduced density matrix is 



|ci(OP 
ct (t)co 



CO 



ci{t)c*o 
+ Efc|cfe(<)p 



In the 



(14) 



which in combination with the condition that |coP + 
Efc|cfe(<)P = 1 - |ci(i)P (Trp(t) = 1) indicates that 
the time dynamics of p{t) is fully determined by ci(t). 

In the limit of T = 0, we can solve for ci{t) exactly 
and obtain the closed form 



ciit)=C- 



ILUoi) 



s[{s + T- iujai)^ + - r2] - Kn'^TTi/T 



where £ ^{F{s)} is the inverse Laplace transform of 
F{s), and K - ia/MiUoi) x ($o/27r)2 « 1/{luqiTi^o)- 
The element pn = |ci(i)p of the density matrix is plot- 
ted in Fig. m for three different values of resistance, as- 
suming that the qubit is in its excited state such that 
Pii(O) = 1. We consider the experimentally relevant limit 
of r ^ ujQi « fi, which corresponds to the weak dissipa- 
tion limit. Since F = l/(2Ci?) the width of the resonance 
in the spectral density shown in Eq. (jlip is smaller for 
larger values of R. Thus, for large R, the RLC environ- 
ment transfers energy resonantly back and forth to the 
qubit and induces Rabi-oscillations with an effective time 
dependent decay rate j{t) = —2^{ci{t)/ci{t)} . 

These environmentally-induced Rabi oscillations are 
a clear signature of the non-Markovian behavior pro- 
duced by the RLC environment, and are completely ab- 
sent in the RL environment because the energy from the 
qubits is quickly dissipated without being temporarily 
stored. These environmentally-induced Rabi oscillations 
are generic features of circuits with resonances in the real 
part of the admittance. The frequency of the Rabi os- 
cillations = sj 7rKfi3/2F is independent of the resis- 
tance since ~ fi\/ ■nL'^C j L\Cq, and has the value of 
^Ra = ^TTfRa « 360 X 10^ rad/scc for Fig.H 
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FIG. 4: Population of the excited state of the qubit as a 
function of time pii{t), with pii(t = 0) = 1 for R = 50 ohms 
(solid curve), 350 ohms (dotted curve), and R = 550 ohms 
(dashed curve), and Li = 3.9nH, L = 2.25pH, C = 2.22pF, 
Co = 4.44pF, Ro = oo and Lq = 0. 



In conclusion, we analyzed decoherence effects in a sin- 
gle superconducting phase qubit coupled to isolation cir- 
cuits with an intrinsic resonance, and emphasized the 
crucial role played by the design of the isolation circuit 
on decoherence properties. In particular, for an RLC iso- 
lation cicuit, we found that the decoherence time of the 
qubit is two orders of magnitude larger than its typical 
low frequency ohmic regime, provided that the frequency 
of the qubit is about two times larger than the resonance 
frequency of the isolation circuit. 

Lastly, we showed that when the qubit frequency is 
close to resonance with the isolation circuit, the non- 
oscillatory Markovian decay of the excited state pop- 
ulation of the qubit, gives in to an oscillatory non- 
Markovian decay, as the phase qubit and its environ- 
ment self-generate Rabi oscillations of characteristic time 
scales shorter than the decoherence time. 
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